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Abstract 

It has recently been shown that by measuring the transverse polarizations of the 
final particles in /i+ — )• e~^e~e~^, it is possible to extract information on the phases 
of the effective couplings leading to this decay. We examine this possibility within 
the context of i2-parity violating Minimal Supersymmetric Standard Model (MSSM) 
in which the /i"*" — >• e~^e~e~^ process can take place at a tree level. We demonstrate 
how a combined analysis of the angular distribution of the emitted electrons and their 
transverse polarization can determine the CP-violating phases of the trilinear R-parity 
violating Yukawa couplings. 
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1 Introduction 



If the neutrino masses are the only sources of Lepton Flavor Violation (LFV), the rates of 
the Lepton Flavor Violating processes such as /i — )■ 67, nN — )■ eA^ and /i — )■ eee will be too 
small to be detectable in the foreseeable future. Any positive signal for such LFV processes 
will be an indication for new physics beyond the Standard Model (SM). Strong experimental 
bounds exist on the rate of these processes ^ and rich literature has been developed on the 
constraints on new physics from these bounds. Currently, the MEG experiment at PSI is 
searching for — )■ 67 and will be able to detect a signal even if Br(|U — )■ 67) is as small as 
0(10-13). 

We can divide the beyond SM scenarios into two classes: (1) Models, such as R-parity 
conserving Minimal Supersymmetric Standard Model (MSSM), within which only even num- 
bers of new particles can appear in each vertex; (2) Models, such as R-parity violating MSSM, 
within which the number of new particles in each vertex can be both even or odd. In case 
of the first class of models, both /i — )■ 67 and /i — )■ eee can take place only at the loop 
level. Within these models Br(|U — )■ eee), being a three body decay, is typically smaller than 
Br(/x — )■ e7) so the latter gives stronger bounds on the LFV parameters (see, however [2]). 
However, within the second class of the models, yU — > eee can take place at tree level and 
its rate can therefore exceed that of /i — )■ e7 which is possible only at a loop level. Within 
the R-parity violating MSSM, it is possible that while Br(yU — e7) is too small to be probed 
even by MEG, Br(/i — )■ eee) is relatively large and close to its present bound [2J. A mild 
improvement on yU — )■ eee can probe R-parity violating MSSM with mass scale well above 10 
TeV which is beyond the reach of the LHC. 

Recently it has been shown that by measuring the transverse polarizations of the final 
particles in yU — )■ e7, /i — e conversion on nuclei and /i — )■ eee, one can measure the CP- 
violating phases of the general effective Lagrangian leading to these processes The 
possibility of deriving information on the CP-violating phases of the R-parity conserving 
MSSM from /i — )■ e7 and the /i — e conversion has been explored in [1]. In the present letter, 
we are going to examine the possibility of deriving the CP-violating phases of the trilinear 
R-parity violating couplings from /i — )■ eee. These CP-violating phases are important as 
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they can be responsible for the creation of the baryon asymmetry in the universe [6]. 

The paper is organized as follows. In section 2, we review the general effective Lagrangian 
that can lead to /i — )■ eee and fi — )■ ej. We outline the information that can be derived 
on the effective couplings without the measurement of the spin of the final particles. In 
section 3, we review the contributions that the effective couplings receive from the i?-parity 
and lepton flavor violating trilinear Yukawa couplings. We then briefly review bounds on 
these couplings. In section 4, we introduce the P-odd asymmetry, A, and the transverse 
polarization. In section 5, we briefly discuss the feasibility of measuring the transverse 
polarization of the final particles. In section 6, we discuss the results that can be derived by 
combined analysis of A and the transverse polarization. Results are summarized in section 7. 



2 Effective Lagrangian in a general framework 

The effective Lagrangian leading to — )■ eee can in general be written as [7] 

£eff = Bi{fxLeR){eReL) + B2ifxReL){eLeR) + 
Ci{jiLeR){eLeR) + C2{jiReL){eReL) + 
GiifJ.RYeRlieR'y^eR) + G2{jiLl''eL){eLlueL) 

-Ar {llLb^^.l^ %eR) {e^e) - Al {Ml^^.lA ^e^) (e7^e) + H.c. (1) 

Notice that by using the identities {<7^)ai3{o'^)js = Se^^e/j^ and {d''^)ai3 = e/35(<7^) 57^70 (where 
Cii = Coo = and eio = — cqi = 1) and employing the fact that the fermions anti-commute, 
we can rewrite the terms on the first line of Eq. ([1]) as 

B B 

— l^{jJ'L'l''eL){eR'-fyeR) - —ifiRYeR){eLlueL) ■ 

This effective Lagrangian leads to [7] 

Brifi^eee) = -L^ + |52p + 8(|Gip + IG2I') 



32^1 



2 mP' TfiP' 

^^ ^[AlGI + ArGI] ^^^ Ji[AlBI + ArBI] 



(2) 
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By studying the energy distribution of the final particles, more information can be derived 
on the effective couplings. For example, let us consider the contributions from and 
which come from a virtual photon exchange. When the invariant mass of an electron positron 
pair goes to zero, the virtual exchanged photon goes on shell. As a result, the corresponding 
Dalitz plot should have a peak at (Pe-+-Pe+)^ = whose height is given by lA^p + IA/jp. The 
combinations that can be derived by studying the energy distributions of the final particles 
are \Gi\^ + + {\Ci\^ + ICsH/lG, \Bi\^ + l^sp, \Al\^ + \Ar\\ RsIAlG; + ArGI] and 
KelAiBl + v4^i?2] (see for example [7j and references therein). By studying the angular 
distributions of the final particles relative to the polarization of the initial muon, one can 
further derive |Gi|2-|G'2p + (-|Ci|2+|C2p)/16, \Bi\^-\B2\\ ReiAiG*-^^^*] 
and Re[Ai5* - ArB;] as well as the CP-odd quantities ImfA^G* + ArGI] and ImfA^E* + 
ARB2]. A combined analysis of angular and energy distribution therefore yields the absolute 
values of all the effective couplings, Bi, B2, A^ and Ar as well as the CP- violating phases 
a.Tg[ALA*j^GlGi], aiglALG^], a.Tg[ALA*j^BlB2] and arg[y4ii?J']. Notice however that there is 
still some information in Eq. ([T]) that cannot be derived by the methods described above. In 
particular, the CP-violating phase arg[i?ii?2] cannot be derived by these methods. Further 
information can be obtained by studying the transverse polarization of e~ in — )• e^e~e~^ 
i. 

Notice that the terms on the last line of Eq. ([T]) come from the effective coupling of the 
photon which gives rise to 

The present bound on this LEV rare decay is very strong Br(/i — )■ 67) < 1.2 x 10"^^ \T\ 
which implies I^lP + \Ar\'^ < 3.3 x 10~^^ GeV^^. Thus, from Eq. (3), we observe that the 
contribution from A^ and Ar to fi eee cannot be larger than 7 x 10"^^ so if Br(/i — eee) 
turns out to be close to its present bound, we can safely neglect the contributions from Al 
and Ar. 
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3 Effects of trilinear R-parity violating couplings 



Within the R-parity conserving MSSM, the slepton mass matrix as well as the trilinear A- 
term involving the sleptons include LFV sources which can induce the effective Lagrangian 
in Eq. ([T]). By relaxing the R-parity conservation, new sources of LFV emerge. In particular, 
the R-parity violating Yukawa couplings in the superpotential 

W = ^LiLjEk with Xijk = -Xjik (4) 

add nine new sources of LFV. That is each nonzero element of Ajj^ is a source of LFV. Some 
of these couplings can contribute to fi —> eee at tree level. Throughout this letter, we will 
focus on the effects of A and set other LFV sources equal to zero for simplicity. The nonzero 
elements of Xijk contain nine phases out of which three can be absorbed by rephasing Lj. In 
this basis, each of the bilinear R-parity terms, ^'^LiHu can be considered as a new source for 
CP-violation. In this letter, we will investigate the possibility of deriving the CP-violating 
phases of Xijk from the transverse polarization of the final particles in — )• eee. 

The contributions to the effective couplings from Xijk have been calculated in [2] and the 
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results follows: 



Al = G2 = Ci 
A321A311 



Co = 0, 



Bo 



2 



(5a) 
(5b) 

(5c) 



a 



127r 



+ 



A321A311 ^ 
A322A312 ^ 



---2 log 



+ 



3m? • 



---21og^ + 



+ 



A323A313 , 



m,- 



2 log 



+ 



a 



247r 



+ 



A321A311 ^ 
A322A312 ^ 



m 
m 



---21og^ + 



2 



m. 



---21og^ + 



A323A313 , 



---21og 



A 



am, 



R 



487r 



^321^311(1 



m 
m 



+ 



3m? ■ 

3m? ■ 
3<^, 



(5d) 



2m?^ 



^) + ^322^312(1 



mt 



2m? 



+ ^323^313(1 



mi 



2<- 



(5e) 



Notice that while the couplings Bi and B2 receive a contribution at a tree level, and Gi 
receive contributions only at the one-loop level. If Xijk are the only sources of LFV, up to 
one-loop level, all other couplings vanish. As discussed in the previous section, the values of 
IE2P, \Al\'^, |v4ij|2, + IC2IVI6 and IG2P + |Ci|Vl6 can be derived by combining 
information from energy and angular distribution of the final particles in /z — )• eee. Thus, 
the predicted pattern for the effective couplings from Xijk can be tested this way. 

Let us now review the various bounds on the couplings from other observations. At 
1-loop level, the Yukawa couplings, Xijk, contribute to the neutrino mass matrix, m,^ [8\ as 
well as to processes such as r — )• fiuu , r — euu and /i — )■ euu P [I] . The upper bounds on 
the components of {m,j)aj3 can be translated into bounds on A's. The prediction of the SM 
for charged lepton decay rates agrees with the measured values so bounds can be set on the 
contribution from A's. In particular, comparing the measured and calculated values of the 
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Coupling(s) 


Bound 


Observable 


Ref 


^133 


94 X 10-^ ( V^^ 




m 


\ 

-^233 


n /I 1 n-4 f "^fR \ 




|1 r\\ 


'^«22 


1-5x10- (,aftv)"'' 




m 


Al2fc 


n no "^^kR 
100 GeV 




m 




n no 

100 GeV 




m 


'^23fc 


n 

"•""^ 100 GeV 




m 


^123-^^32 


o in-7 ™Ah '"fH 
* ^ (,(100 GeV)^ ) 




m 



Table 1: Bounds on the trilinear i?-parity violating couplings. The masses of e^, hr and tr 
are respectively indicated by mg^^, rrie^^ and rrie^^. 



ratio Rr = T{t — )■ et'e'^r)/r(T — ^ fJ'J^fj.^r) gives bounds on A's. In our analysis, we pick up 
values of A that respect these bounds. The bounds that we use are summarized in Table 1. 
The third column shows the observable from which the bound is extracted. Notice that the 
bound on Ai2fc comes from the numerical value of the CKM matrix. At first sight, this might 
seem counterintuitive. Remember however that the value of Vud is extracted by comparing 
r(n -T- peue) with r(/i — )■ ez/gt'^). The unitarity condition of the CKM matrix combined by 
the extracted values of the CKM matrix elements yields a bound on the contribution from 
^2ik- For a full review of the bounds see [TOj . 



4 Transverse polarization and asymmetry 

Let us define P-odd asymmetry. A, as 

Jq {dV /dcos9)dcos9 — J^^{dT /d cos 9) d cos 9 

A = 1 . (6) 

J_^{dT/d cos 9)d cos 9 

A nonzero A violates parity. In fact, A is sensitive to the P-odd combinations such as 
|-Bi|^ — |-B2p- The polarizations of the final electron in /x"*" — )■ e^e~e^ can be defined as 

dr I dr \ 



+ , 5*4- dcos0ls„--T. = s„_-T. /-^ 



(^t- ) - — — (7) 



,s 4- ,s dcos 8 



where i = 1,2,3. T3" is the longitudinal direction, T3" = Pe-/|-Pe-|- ^i" ^^'^ ^^i^ 
vectors in the transverse directions defined as 

f, ^ and f, ^ (8) 

\s^y<Pe-\ \T2XP,-\ 

Finally, 6 is the angle between the momentum of e~ and the polarization of the muon: 

6 = arccos[s^ ■ Pe~/{\s^\ ■ \Pe-\)] ■ 

As shown in |3j, while (sf-) is sensitive only to the absolute values of the couplings, 
the transverse polarizations (sf-) and (sf-) are sensitive to the CP- violating phases in the 
effective Lagrangian. Straightforward but cumbersome calculation shows that within the 
present model with the coupling pattern in Eqs. ([5]), 

{sj,-) oc ^BiB; - UAtiBl] 

and 

{st-) ^ Q[BiB; - UAnBl] . 

Of course, similar formulas hold for the transverse polarization of the final positron in /i~ — )■ 
e~e'^e~. 

Let us suppose Br(/i — )■ eee) is measured and found to be close to 10~^^. Let us moreover 
suppose that the Dalitz plots reveal that Bi and B2 give the dominant contribution to 
this decay as it is expected in the case that LFV effects originate from Xijk- If the MEG 
experiment at PSI reports a null result {i.e., Br(/i — > 67) < 10^^'^ ), within the context of 
MSSM (to be tested at the LHC), such a set of conditions attests our assumption that Xijk 
is the prime source for LFV. If MEG also finds a signal for fi — > ej, other LFV terms, such 
as slepton masses or the trilinear soft supersymmetry breaking terms, might contribute to 
Br(/x — 67) but as we discussed earlier, even in this case, we can neglect the contribution from 
Al and Aji to Br(/i — )■ eee). LFV terms in slepton masses or the trilinear soft supersymmetry 
breaking terms can contribute to other effective couplings in ([T]) but the effect will be loop 
suppressed and negligible. Notice that although within the model under our study, Gi and 
Aji are given by the same combinations of Xijk, Gi is enhanced by a factor of log(m?/m^) so 
we do not neglect its contribution [see Eq. (^]. Neglecting the effects of A^, the formulas 
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for A, (sf-) and (sf-) will have forms: 

3(|Bi|2+ |52P + 8|Gi|2 



and 



/3 „ ^ 4g?[Eii?*]P.sing 

^ ^i"^ |5i|2(3-P^cos^) + |52|2(3 + P^cose)+24|Gi|2(l-P^cos^)' ^ ^ 

|fii|2(3-P^cose) + |52P(3 + P^cose)+24|Gi|2(l-P^cos^) ' ^ ^ 

in which P^ is the polarization of the initial muons. 
Moreover the longitudinal polarization is given by 

Notice when the electron is emitted in the direction perpendicular to the spin of the 
muon, the transverse polarization is maximal. In our analysis, we will set Q = 7r/2 which 
experimentally means we study the data from the polarimeter collecting electrons with mo- 
mentum perpendicular to s^. If more than a single polarimeter is installed, more data can 
of course be collected. The average polarization can be defined as 

I-ii^f^) [dT {jj, ^ eee) / d cos 9] dcosO 
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{sf-) = \ ' — ^ . (li: 

/_-^[(ir(yU — )■ eee)/(icos^] (icos^ 



For i = 1,2, (sf-) = f P/x(sj.-)|6»=|- It is noteworthy that if the values of Br(/i — eee), A 
and {sf-)\0=2L are measured, the numerical values of |i?i|, |i?2| and \Gi\ can be derived. 
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5 Feasibility of measuring the transverse polarization 

The typical experimental setups devoted to the study of muon decay (such as the MEG 
experiment or the experiment at TRIUMF described in [13]) can be summarized as follows. 
A proton beam collides on a target producing pions. Charged pions are stopped in the target 
and decay at rest into a neutrino and a muon. Muons at production are 100 % polarized up 
to neghgible correction due to the neutrino mass Muons exit the first target and are 
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transmitted to a second target where they stop and then decay at rest. Based on the setup of 
the experiment, muons of either positive or negative sign can be selected to be transmitted 
to the second target. Negative muons would form bound states with atoms in the second 
target so we focus on the decay of positive muons which decay as free states. When muons 
decay, they are still highly polarized. The degree of depolarization from the production to 
decay depends on the setup of experiment. Especially if only the muons produced at the 
surface of the first target are collected (as it is done both in MEG [12j and in the experiment 
described in pJij), the depolarization will be quite negligible. For example, in the TRIUMF 
experiment described in [13], the polarization remains above 99% until the muon decay in 
the second target. For the purpose of the present analysis we can safely replace = 100%. 
Notice that unlike the case of [T3], even a moderate accuracy in knowledge of will be 
sufficient to perform the present analysis so from this aspect, it is easier to carry out the 
present measurement |13j . 

Let us now discuss the possibility of measuring the transverse polarization of the electron 
in /i+ — )■ e^e~e^ . First, let us recall the measurement of the transverse polarization of the 
positron in /i"*" — i- e^VeV^ which was carried out about 25 years ago to extract the Michel 
parameters P^. The outline of the polarization measurement is as follows. The emitted 
positrons interact with electrons in a target which are polarized in a direction perpendicular 
to the momentum of the emitted positrons by a magnetic field. The electron positron pair 
annihilate into a photon pair. By studying the azimuthal distribution of the final photon 
pairs, the transverse polarization of the emitted positron can be determined. In our case, 
we are interested in measuring the transverse polarization of the final electron rather than 
the final positron. Fortunately, a similar setup can be employed to determine the transverse 
polarization of the emitted electron. Of course in this case instead of annihilation into a 
photon pair, the Moller scattering (e]"e^ — )■ ejej) will take place. Similar to the case of 
e~e^ — )■ 77, the azimuthal distribution of the final particles carry information on the spin 
of the emitted electron. Let us take the z direction to be parallel to the momentum of the 
emitted electron (the momentum of e]~) and take x to be in the direction of the spin of the 
electron at rest in the magnetized target {i.e., the spin of ): 

Pi = (^A;2 + m2,0,0,A;) and P2 = ("^e, 0, 0, 0) . 
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The spin of ei can be described by (a b) where |ap + |6p = 1: 

(s^) = 2^[ba*], {sy) = 2^ba*] and {s,) = |a|2 - |6p . 
The angular distribution of the final particles are described by {0,(f)): 

P3 = + ml, k' sin 9 cos 0, k' sin 6* sin 0, k' cos 6*) and P4 = Pi + P2 — -fs , 

where the energy-momentum conservation implies 

, , rupk 
k' 



me + k{l — cos 6) 
A cumbersome but straightforward calculation shows that 



2- da 



dcos9 d(f) 



cos (f)d(f)— (12) 



2 I |2xA;'V(l -cos^)sin^ -cos^) - cos ^) - me(l - 2A;' cos ^/A;) \ 

' ' ~ ' 647rmf (A; - k') \ {k - k'f ^ {k - m,f {k - k'){k - m^) J 
and 

f27r 

— - sin 2(/)d(/) = (13) 

iQ dcosO ^'^ 

k'^e'^{l — cos9) sin^ 9 f rrie/k rrie/k 1 



^0 



MTimiik-k') \{k-k'f {k ~ m^)'^ {k - k'){k - rrie) _ 
Thus, by measuring the partial cross section, da/ {d cos 9 d(f)), and taking the above integrals, 
a and b (up to an overall phase) can be determined and the spin of ei can be therefore 
reconstructed. 

There are two problems that complicate the measurement: (1) In the lab frame where ej 
is at rest, the majority of the electrons are scattered in the forward direction within a narrow 
cone with opening angle of 0{2me/k) where k ~ ~ 30 MeV. The same problem existed 

in the case of measuring the polarization of the emitted positron. (2) The scattered electron 
63 can again scatter on the electrons in the magnetized target before exiting it. Multiple 
scattering will distort the azimuthal distribution in which the information of the spin of the 
initial electron is imprinted. The same problem existed is the case of measuring the spin of 
the positron as the photons produced in the electron positron annihilation could Compton 
scatter on the electrons in the magnetized target. In both cases, the total scattering cross 
section is of order of e^/ {IGnrriek). Fortunately, there are established techniques to overcome 
these difficulties. 
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6 Combined analysis of A and {sf-) 

By rephasing the leptonic fields, three out of nine phases in the A couphngs can be absorbed. 
Let us consider the basis in which A311, A211 and A312 are all real. The rest of A's in this 
basis can in general be complex among which the phases of A321, A212, A322 and A323 can lead 
to a nonzero (sf-). We first concentrate on the phases of A321 and A212 which contribute to 
Bi at a tree level. We then comment on the rest of phases. As we discussed in the previous 
section, we can safely replace = 100 % and that is what we do in the following. For any 
given polarization, our results can be simply rescaled. 

Figs. 1 and 2 demonstrate the dependence of (sf-) and A on the phases. To draw these 
scatter plots, we have assigned random numbers at a logarithmic scale from 10~^ to the 
upper bound on |A|'s. (|A|'s take up random values in a logarithmic scale.) As explained 
in the caption, various values are assigned to the phase of either A321 or A212, setting the 
other one (as well as the rest of phases) equal to zero. The phase of Xijk is denoted by 
(Pijk- We have selected the configurations of A for which Br(/x eee) lies within the range 
Br(/x — >■ eee) = 5 x 10~^^(1 ± 10%). From an experimental perspective, this means that 
we have supposed Br(// — > eee) is measured to be in the range 5 x 10"^^(1 ± 10%). The 
10% uncertainty is a nominal value that we have taken as an example to highlight the fact 
that the branching ratio measurement will suffer from a finite uncertainty. Our results are 
robust against varying the value of this uncertainty. In fact by rescaling the coupling by a 
5N percent, Br(^ — )■ eee) will change by 4:SN% but (-Sj.-) and A, being ratios, will not vary. 

Points denoted by pink crosses are the points at which (^321 — and (^212 takes up random 
values in (0, 2tt). On the contrary, those shown by green circles correspond to the cases that 
^212 — and (fi32i takes random values in the range (0, 27r). As seen from the figures, the 
areas over which pink x and green o are scattered completely overlap which means if |A|'s 
are unknown, the two solutions are indistinguishable. However, valuable information from 
{sf- ) can be derived. For example if both CP- violating phases vanish, {sf- ) also vanishes 
(points indicated by green V). 

Remember that, while Bi and B2 receive nonzero contributions at a tree level, Gi receives 
a contribution only at a loop level. As long as Gi has a neghgible value, for given A and 



12 



Br(/x — )• eee), \Bi\ and I-B2I are fixed. It is straightforward to check that, for — )■ 0, 

\A\<^ and \{sf-)\ < . (14) 

As seen in the figures, the majority of points lie inside an oval-shaped region whose bound- 
aries are given by Eq. (|T^ . These are the points for which the contribution from loop 
suppressed can be neglected. As seen from Fig. 1, there are only few points (about 

2 percent of all points) lying outside the oval-shaped region. At points with A < —1/3, 
the contribution from dominates. For A < —0.5, we find Ks^^-)], ^ 0.1 which 

is expected from Eqs. (jS]) and (|TOl) . In Fig. 2, we have removed the points for which 
SlGip/d-BiP + I-B2P) > 0.05. As a result. Fig. 2 does not include points outside the oval- 
shaped region. If the pair {A, {sf-)) turns out to be outside the oval-shaped region, it means 
IGil is relatively large. This can happen if A322 is more than 50 times larger than the rest 
of Xijk (see Eqs. and table 1) so ^ < —1/3 indicates that the flavor structure of the Xijk 
coupling should be hierarchical. 

As mentioned above, in the limit Gi — ?■ 0, for a given A and Br{fi — eee), \Bi\ and 
|i?2| are fixed so {sf-) determines arg[i?2-Bi]. Thus, if all the phases except (^9321 are zero, 
(s^-) will determine v^32i. This can be seen in Fig. 2. That is impressive that ^9321 can 
be derived even without knowledge of |Ajjfc|'s (of course under the assumption of a single 
nonzero phase). Deriving the value of ^9212 is going to be more challenging even when we set 
all the other phases equal to zero. As seen from Figs. 1 and 2, for a given A and a certain 
value of (p2i2, depending on the configuration of the absolute values of the A's, {sf~) can 
take any value between zero and its maximal value which is (2/3) (1 — 9 AY^'^ sin Lp2i2- This 
is understandable as B2 receives contributions from various terms (compare Eqs. fl5b|) and 
( 15c|) ) so unlike the case of (^5212 = 0, ip32i 7^ 0), in this case, the nonzero phase is not given 
byarg[5i5*]. 

From Fig. 2, we observe that for |^| < 0.2, by simultaneous measurements of A and 
{sf-) with reasonable accuracy, even without independent knowledge of \Xijk\, solutions with 
(v'321 = 7r/2,y92i2 = 0) and {(p32i = 7r/4,(y92i2 = 0) can be distinguished (see points denoted 
by violet > and pink □). Notice that all points denoted by red A and blue < corresponding 
respectively to (V5321 = 0,^9212 = 7r/2) and {(p32i = 0,^9212 = 7r/4) lie above the horizontal 
axis. Solutions with positive and negative sin (^9212 are distinguishable but deriving the value 
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of v?2i2 without knowledge of \Xijk\ does not seem to be practical. We have found that the 
contributions from the phases that enter only at loop level {i.e., (p322, V'sis and V5323) to {sf~) 
are smaller than 0.1. Thus, in deriving the values of 1^212 and ^9321 from (sf-), the potential 
contributions from the rest of the phases can be ignored. 

As discussed above, an upper bound on IGil can considerably simplify the analysis and 
solve the degeneracies. Although jGip (more precisely, 1(711^ + IC2P/I6) can in principle be 
extracted from the energy distribution of final particles, within the present model, its value 
will most probably be too small to be measured so in practice only an upper bound on 
can be extracted as we have assumed in deriving Fig. 2. Notice that 

IGil^O A 

That is while if Gi and G2 (or Ci and C2) gave the main contribution to — )■ eee, we would 
expect that {sf-)\e=^ /A = 1. The ratio of longitudinal polarization to A can therefore be 
regarded as a cross-check for the smallness of |Gi|. 

To draw Figs 1 and 2, we have used the spectrum at the a benchmark [TJ] as the input: 
i.e., We have set m^^ = m^^ = 285 GeV. For two reasons, we expect the results to be 
robust against varying the input masses: (i) Varying m?^ and mf,^ is respectively equivalent 
to rescaling A211 and A311 (see Eqs. (l5b|) and (|5cl)). (ii) Both A and (sf-) are defined as 
ratios so the dependence on the supersymmetry scale disappears. We have re-drawn the 
diagram for different benchmarks. As expected, the results are not sensitive to the input for 
the mass spectrum. 

It is noteworthy that if the only sources of LFV are the Xtj^s giving rise to Br(yU — )■ eee), 
Br(/i — )■ e7) will be smaller than 10^^^ so if the MEG experiment reports a /i — )■ e7 signal, 
sources for — e7 other than Ajj^'s must exist. 

As seen in Figs. 1 and 2, both A and (sf-) can be relatively large so as long as the 
errors in their measurement {i.e., 6A and 6{sf~)) are below ~ 0.1, their nonzero values 
can be established. Suppose the numbers of electrons studied to derive A and {sf~) are 
respectively A/a and Ms. We roughly expect the statistical errors to be 6 A ~ l/^/J\^A 
and S{sf-) ~ l/vCV^. To derive A, the majority of the emitted electrons can in principle 
be employed, so with a few hundred /x+ — )■ e+e'e"*" decays, the statistical error in the 
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measurement of A will be reasonably small and below 0.1. However, we expect only a 
fraction of the emitted electrons, r, to enter the polarimeters so for establishing nonzero 
(sf-), the total number of /i+ e^e^e^ decays has to be larger than 100/r. That is if 
r ~ 10%, the total number of ii^ — )■ e^e~e^ has to be larger than a few thousand. 

In the above analysis, we have employed information on the i?-parity violating couplings 
from only the LFV rare decays. The R-parity violating couplings can in principle be directly 
measured by accelerators. |Ajii| leads to a resonant production of z>j in a e~e^ collider 
(e~e"'" — )■ Ui) so |Ajii| can be derived provided that |Aiii| > 10^^ and the center of mass 
energy is equal to the mass of z/j [161 HO]- Moreover, the A couplings can lead to Vi — t- l^li , 
Xi — ^il^lk ^"^^ xt ~^ ^k^il^ 1 ^t^k^J [13 where Xi and xt a-^e the lightest neutralino 
and chargino. Thus, by measuring the decay length and the flavor of the final charged 
leptons, \\ijk\ can in principle be extracted. If \\ijk\ < 0(10^^), the decay length will be 
too small to be resolved [10]. This method is therefore sensitive only to the values of \\ijk\ 
smaller than (9(10~^)1^ If |A|'s are all smaller than 10~^, each \\ijk\ might be extracted 
from Vi — > Xi ~^ ^i^~jK ^"^^ xt ~^ ^ki^i^j' Jt^t^J ^^^^ '^^^^ Br{fi — )■ eee) will 

be too small {Br{fi — )■ eee) < 10^^^). Let us now consider the range, nisusy ~ 100 GeV, 
A2115A311 > 10^'^ and Xijk ~ 10~^ with ijk 7^ 211,311. In this range, e'e"*" — z>j yields 
|Ajii| and Br(yU — > eee) is close to the present bound. Moreover, for ijk 7^ 211,311, Ui — )■ 
l^lk and Xi ^i^~jK ^^^^ have a resolvable decay length but the point is that the decay 
modes involving A211 and A311 will dominate and lead to a decay length too short to be 
observable: e.g., — )■ e+e~)/r(z>^ — )■ r+r^) ~ 10^. Thus, even in case that the i?-parity 
conserving decay modes are kinematically forbidden, extracting the decay lengths will be 
quite challenging. Let us however suppose that these experimental difficulties are partly 
solved and certain \\ijk\ (but not necessarily all) are measured. Such achievement might 
not be out of reach if Aju ~ 10~^ and the rest of A's are of order of 10~^. Complementary 
information can then be derived from Br(/i — )■ eee), A and (srp-): Neglecting the effects. 



^Even if the decay length is not resolved, a combination of |A|'s may be extracted. For example, consider 
chain processes e~e~^ — )■ Z* — !■ i>ii>i and the subsequent decays i>i — > e~ fi~^ and Vi — > fi~ fj,'^. Such a chain, 
being LFV, is not contaminated by the SM or i?-parity conserving MSSM so even if the decay lengths of 
i>i — )■ fj,'^ and i>i fJ- fJ-^ sue too short to be resolved, the possibility of deriving information on the 
relevant couplings is still open. Considering all such possibilities is beyond the scope of the present letter. 
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Br(/x — )• eee) and A give |-Bi| and I-B2I which to leading order correspond to IA311I ■ IA321I and 
1^211^212 + i^l /ttT'I )AqiiA312|, respectively. Thus, if A311 is extracted from e+e~ — )■ at 
ILC, the measurement of Br(/i — > eee) and A gives IA321I. If IA211I, IA212I, IA311I and IA312I 
are all derived by accelerators, this method will give the phase of A2i2- The measurement of 
(sf-) will then yield the phase of A321. 

7 Conclusions and discussion 

The trilinear i?-parity violating Yukawa couplings, XijkLiLjEkf^ can lead to /i — )■ eee. In 
particular, A321, A311, A211, A312 and A212 contribute to /i — )■ eee at tree level. By rephasing 
the lepton fields, we can go to a basis in which A311, A312 and A211 are real. This exhausts 
the freedom to rephase the other fields so A321 and A212 can in general be complex and can 
be considered as sources for CP-violation. Their phases can induce transverse polarization 
for e~ in the direction T2 = x Pe-/|s)j x Pe-\. Thus, by measuring this polarization, 
one can derive information on the CP-violating phases. We have shown that for maximal 
CP-violation, can reach 2/3 so with a moderate sensitivity to (sf-), CP- violation 

can be established. The sign of the CP-violating phase can also be determined by measuring 
the sign of {sf-). 

We have also studied the P-odd asymmetry, A defined in Eq. and discussed the 
information that from a combined analysis of A and (sf-) can be obtained. For the majority 
of the A configurations, the tree level effects dominate so the effective coupling is much 
smaller than \Bi\ and IP2I and its effects can therefore be neglected. In this case, \Gi\ 
will be too small to be measured but an upper bound can be put on by studying 
the energy distribution of the final particles in /i — )■ eee or as discussed in the present 
paper by combining information on A, Bi{fi — )■ eee) and (sf-). We have noticed that if 
8|Gip < 0.05(|-Bip + I-B2P), the analysis becomes much simpler because of two reasons: (i) 
The contributions of the phases of A322, A323 and A313 to (sf-) can be neglected. These are 
the couplings that contribute only at a loop level, (ii) For given A and Br(/i — )■ eee), the 
absolute values of Bi and B2 are fixed which simplifies the analysis. In particular, restricting 
the analysis to a single CP- violating phase, the simultaneous measurement of A and (sf-) 
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yields the phase of A321 even if the values of \Xijk\ are not a priori known. Degeneracies 
however exist between solutions for which the phases of A212 and A321 are both nonzero. By 
measuring A and {sf-), different classes of solutions can be distinguished. If the absolute 
values of A are measured by an accelerator-based experiment (or by some other methods), 
the measurements in /i — )■ eee can yield both phases. In fact, if the tree level contribution 
dominates, the relevant A parameters can be over-constrained. 

Neglecting the contribution of |Gi|, we find -1/3 < ^ < 1/3 and -2a/(1 - 9A'^)/9 < 
{sf^) < 2a/ (1 — 9A'^)/9. Only at a small fraction of the parameter space where (A322 ^ 
rest of A), the loop effects can dominate and lead to ^ < —1/3. Within the model un- 
der consideration, A < —1/3 therefore indicates a hierarchical flavor structure for the A 
parameters. 

Stopped yU~ would form bound states with the atoms before they decay. Because of this 
technical difficulty, we have concentrated on the decay of fi^ and the spin of the flnal electron 
emitted from it. Similar consideration holds for the positron in free decay of negative muon; 
i.e., for e"*" in /i~ — )■ e~e~^e~. Within the present model, the transverse polarizations of the 
electrons in /i~ — )■ e^e+e^ (or that of the positrons in /i"*" — > e^e~e"^) are loop-suppressed. 

There are established techniques to measure the transverse polarization of the positron 
based on the azimuthal distribution of the photon pair produced by the annihilation of the 
positron on the polarized electrons in a thin magnetized target [H]. We have shown that a 
similar setup can be employed to measure the transverse polarization of the electrons, too. In 
fact, the azimuthal distribution of the flnal electrons in MoUer scattering, e^e^ — )■ Cg ej with 
polarized 63 is sensitive to the polarization of e]~. The challenges before this measurement 
are similar to the ones in [H] and can be overcome by similar methods. 

One can repeat similar discussion for three body LFV decays of r lepton such as r — )■ fififi 
or r — )■ eee. The measurements of the angular distribution and polarization of the flnal 
particles in these decays can provide complementary information on the Xijk couplings. Such 
a study will be presented elsewhere. 
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Figure 1: Transverse polarization of the electron in yU"*" — t- e"'"e~e"'" versus the P-odd asym- 
metry A defined in Eq. ([6]). The input values for masses correspond to the ones at the a 
benchmark in [15] with mj>^ ~ rrij^^ ^ 285 GeV. Random values between 10~^ up to bounds in 
Table 1 are assigned to the A couplings and points at which Br(yU — )■ eee) G 5- 10^^^(1±10%) 
are selected. To calculate {sf^) and A, we have set P^j = 100% and 6* = 7r/2 (see definitions 
in Eqs. flTfS]) ). ipijk is the phase of Xijk- Points with different colors and symbols correspond 
to a nonzero value for ipijk as described in the legend. For each set, the rest of phases vanish. 
For points shown by x and o, the corresponding phase takes random values at a linear scale 
from (0,27r). 
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Figure 2: The same as Fig.lH except that we have removed the points at which 8|G'ip/(|i?ip+ 
iSaP) > 0.05. 
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